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FLUTTER, VIBRATION, AND BUCKLING 
OF TRUNCATED ORTHOTROPIC CONICAL SHELLS WITH 
GENERALIZED ELASTIC EDGE RESTRAINT* 

By Sidney C. Dixon and M. Latrelle Hudson 
Langley Research Center 

SUMMARY 

A theoretical investigation has been made of the flutter, vibration, and buckling char- 
acteristics of orthotropic truncated conical shells with generalized elastic edge restraint. 
The problem is simplified by making certain assumptions which are considered justified 
by the results for cylinders; the shell analysis is of the classical Donnell type, in-plane 
inertias and structural damping are neglected, and the aerodynamic loading is represented 
by the inviscid two-dimensional quasi-steady approximation (modified piston theory). The 
principle of virtual work is utilized to formulate the problem. An approximate solution is 
obtained by the generalized Galerkin method in terms of assumed displacement functions 
that must satisfy only geometric constraints. A stress function is introduced so that only 
the normal displacement need be assumed. The equation for the stress function is solved 
exactly in terms of the arbitrary coefficients in the assumed displacement. 

The accuracy and limitations of the present analysis are illustrated by comparing 
numerical results for buckling and vibration with results of other investigations for vari- 
ous boundary conditions, applied loads, and shell geometries and stiffnesses. Numerical 
results for flutter of conical shells reveal that trends established for buckling and vibra- 
tion do not necessarily hold for flutter. Increasing cone semivertex angle (up to about 
15°) for length-radius ratios greater than about 1 increases the resistance to flutter even 
though the buckling load and natural frequencies decrease. Also, increasing the axial edge 
restraint can decrease the resistance to flutter even though both the buckling load and 
natural frequencies are increased. These phenomena are attributed to the fact that flutter 
depends on frequency spectrum as well as magnitude. Sufficient numerical results are 
presented to permit the determination of the flutter condition for simply supported iso- 
tropic conical shells for a wide range of cone angle, length-radius ratio, and radius- 
thickness ratio. Results are also presented to indicate some effects of variations in axial, 
radial, circumferential, and rotational edge restraint, applied load, and ring or stringer 
stiffeners. 

*Part of the information presented herein was included in a thesis entitled 
"Buckling of Orthotropic Truncated Conical Shells With Elastic Edge Restraint Subjected 
to Lateral Pressure and Axial Load" submitted by Sidney C. Dixon in partial fulfillment 
of the requirements for the Degree of Doctor of Philosophy in Engineering Mechanics, 
Virginia Polytechnic Institute, Blacksburg, Virginia, March 1967. 



INTRODUCTION 


Prior to about 1950 there was a paucity of theoretical and experimental informa- 
tion on the static and dynamic behavior of truncated conical shells primarily because of 
the mathematical difficulties encountered in theoretical analyses and because thin conical 
shells had limited application as structural components. However, more recently thin 
conical shells have been used as load-carrying structural components in adapter sections 
of multistage rockets and in proposed truncated-cone nose-cap configurations for plane- 
tary entry vehicles; therefore, a thorough knowledge of the static and dynamic character- 
istics of such configurations is required. This need and the availability of digital com- 
puters to aid in circumventing the mathematical difficulties encountered in cone analyses 
have resulted in considerable research effort on the buckling characteristics of cones 
(for example, refs. 1 to 4) and on the vibrational characteristics of cones (for example, 
refs. 5 to 7); however, only limited results are available for the effects of various classi- 
cal boundary conditions or elastic edge support, which have been shown to be important 
for cylindrical shells (for example, refs. 8 to 13). Also, there is still an almost complete 
dearth of theoretical and experimental information on flutter of conical shells (ref. 14). 

The present investigation considers the flutter, vibration, and buckling of a trun- 
cated orthotropic conical shell with generalized elastic edge restraint. The problem is 
simplified by making certain assumptions which are considered justified by the results 
for flutter of cylinders. The principle of virtual work is utilized to formulate the gov- 
erning equations of the problem. An approximate solution is obtained by the generalized 
Galerkin method (refs. 15 and 16); that is, a procedure in which the displacement func- 
tions must satisfy only geometric constraints. A stress function is introduced so that 
only the normal displacement need be assumed. The equation for the stress function is 
solved exactly in terms of the arbitrary coefficients in the assumed displacement. 

The accuracy and limitations of the present analysis are illustrated in the appen- 
dixes by comparing numerical results for buckling and vibration with results of other 
investigations for various boundary conditions, applied loads, and shell geometries and 
stiffnesses. Flutter results are presented to indicate the effects of conicity and shell 
geometry on the flutter characteristics of simply supported isotropic conical shells. 
Sufficient numerical results are presented in the form of general flutter boundaries to 
permit the determination of the flutter condition for wide ranges of cone semivertex angle, 
length-radius ratio, and radius-thickness ratio. Results are presented to indicate the 
effects of applied loads, ring and stringer stiffeners, and finite elastic edge restraint 
provided by massless uncoupled springs; the limiting cases of free and clamped edges are 
also considered. 
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Certain lengthy expressions developed in this analysis are required for complete 
documentation. Since they are of limited interest, these expressions are included in a 
"Supplement to NASA TN D-5759,” which is available upon request. A request form is 
included at the back of this paper. 


SYMBOLS 


The units for the physical quantities defined in this paper are given both in the 
U.S. Customary Units and in the International System of Units, SI. Appendix A presents 
factors relating these two systems of units. 


A 

a o’ a o> a m 

B 

By, Be 

Bye 


cross-sectional area of stiffener 

coefficients in assumed normal deflection function 

extensional stiffness, — — — 

1 - M 2 

extensional stiffness in y- and e-direction, respectively 
in-plane shear stiffness 



c air speed of sound in air 

c m speed of sound in cone material 

'F'h 2 

D bending stiffness, — r 

12(1 - jU 2 ) 

Dy,Dg bending stiffness in y- and e-direction, respectively 

Dyg in-plane twisting stiffness 

d stiffener spacing 
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Ill I 


II I 


I I I 


l I 


E Young's modulus 

F stress function 

F U ,F V ,F W forces in direction of u, v, and w applied to cone edge by elastic medium 
f = e -x F 

G shear modulus 



I moment of inertia of stiffener about its centroid 


i,j integers 

J torsional constant for stiffener 


Kij stiffness coefficient representing restraint of elastic medium at cone edges 

stiffness of springs resisting axial, circumferential, and radial edge dis- 
placements, and edge rotation, respectively 

ki,k 3 stiffness of springs resisting meridional and normal edge displacement, 

respectively 

L height of conical frustum 

L x ( ),L 2 ( ),L 3 ( ) operators defined by equations (B15), (B33), and (B59) 
l slant length of conical frustum 

M maximum number of terms in sine series 
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local Mach number 


M t 


Moq 

My,M 0 ,My fl 

My,Myg 

M w 

m ,p 

Ny ,Ng,Ny Q 

Ny ?Ny Q 

n = — ^ — 
sin a 


free-stream Mach number 
moment resultants 

externally applied moment resultants 

moment applied to shell edge by elastic medium 

integers 

force resultants 

externally applied force resultants 


P applied axial load, positive for compression 

Pj,P 2 constants defined by equation (B53) 

p lateral pressure load on shell (positive for external pressure) 

p g surface loading induced by shell vibration or flutter 

p^. total surface load acting on shell, P + P s 

Qy ,Qq transverse shear force resultants 
Qy externally applied transverse shear force resultant 

q local dynamic pressure of airflow 

R radius of conical frustum normal to axis of revolution 


radius of cross section of circular end ring 
v3 Kq 

V 


S il=(c^3-) 3 § 0.i=1.2,3) 
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S 4r(c3Ts) ^ ( S 4j= S j4. J- 1,2.3) 

Ri K44 
S 44 = cos a d^" 

t time 

U s strain energy of shell 

u,v,w displacements of middle surface of conical frustum (see fig. 1) 

U H’ U V radial and axial components of displacements 

x transformed coordinate, x = In — 

’ y 2 

y l 

x, = In — 

1 y 2 

y, z orthogonal coordinates (see fig. 1) 

z distance from centroid of stiffener to middle surface of conical shell, 

z c + z s 

z c distance from neutral axis of sheet-stiffener combination to middle surface 

of conical shell 

z s distance from neutral axis of sheet-stiffener combination to centroid 

of stiffener 


a semivertex angle of cone 



y 

A 

50 e 


shell mass per unit area 
frequency spectrum parameter, 



virtual work done by generalized elastic edge restraint 


6S2j virtual work of boundary forces 
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6£2p virtual work of total surface pressure 

total virtual work of all external forces 
strains at middle surface 
V number of circumferential waves 
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X 

X = 


circumferential coordinate (see fig. 1) 


flutter parameter, 


2qRj 


(l - M 2 )(Ri/h) J 


where x = tanh/ 0.5 ^ 


V R i, 


'-o - 


x l> x 2 




B, 


constants defined by equations (B55) 


Poisson's ratio for isotropic cone 




MQjM' 


Poisson's ratio for bending in 6- and y-direction, respectively 
y Poisson's ratio for extension in 6- and y-direction, respectively 
p local density of air 

r thickness of end ring 


< p = 6 sin a 

X0>Xy>Xyg curvatures at middle surface 
circular frequency 


u> 


Ct) r =| 


Dy\V2 

A 


U)i 


first natural frequency for given value of n 



second natural frequency for given value of n 


“2 

Superscript: 

* indicates variable is function of x only 

Subscripts: 

A conditions existing prior to buckling, vibration, or flutter 

B conditions existing at buckling, vibration, or flutter 

cr critical value required for buckling 

cyl equivalent cylinder 

0 conditions at 77 = 0 

Im imaginary part 

Re real part 

R ring 

S stringer 

1 small end of cone 

2 large end of cone 

A subscript preceded by a comma indicates partial differentiation with respect to 
that subscript. 


ANALYSIS 

Statement of Problem 

The flutter, vibration, and buckling characteristics of orthotropic truncated conical 
shells with generalized elastic edge restraint are investigated. The configuration and 
the notation for coordinates and displacements, force and moment resultants, and loads 
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are shown in figure 1. The cone is subjected to uniform supersonic airflow over its 
external surface at a local Mach number Mj, dynamic pressure q, and density p. 

Also, the cone is subjected to either internal or external lateral (normal) pressure 
loading p and an applied axial load P. The pressure is uniform and acts in a constant 
direction. The movement of the cone edges is resisted by a general elastic medium that 
can be represented by a 4 by 4 matrix of stiffness coefficients Kjj (ref. 17). The mass 
effects of the edge support are represented by a 4 by 4 matrix of inertia coefficients Gij. 
The off-diagonal terms in these boundary-condition matrices account for the coupling that 
can exist between the various displacements due to the interaction of the shell and the 
elastic medium. This coupling must be considered if the edge restraint is provided by a 
ring (ref. 18). 

The shell is stiffened in such a manner that the cone material has orthotropic stiff- 
ness properties in two orthogonal directions. The principal axes of orthotropy are alined 
with the shell coordinate axes (0, y). The five extensional stiffness properties By, Bg, 
By q, p^., and p'g and the five bending stiffness properties Dy, D$, Dyg, Py, and 
Pg are needed to describe the material properties. From the reciprocal theorem (for 
example, ref. 19), the following relations must hold: 

Bylx'g = 

( w 
DyP# - DgPyl 

Thus, there are eight independent elastic constants. 

The flutter, vibration, and buckling characteristics are determined for various com- 
binations of cone geometry, stiffness properties, elastic edge restraint, and loadings p 
and P. 


Assumptions 

The stability characteristics of truncated conical shells should be basically similar 
to the stability characteristics of cylindrical shells. Therefore, the wealth of information 
on cylinder buckling and flutter should be a useful guide as to what simplifying assump- 
tions can be made without seriously affecting the accuracy of the theoretical results. 

Only flutter resulting from coalescence of the first two longitudinal modes of 
vibration is considered since it has been shown theoretically (ref. 20) and experi- 

mentally (ref. 21) that this type of flutter represents the critical stability boundary for 
cylinders. This type of flutter is usually associated with a mode shape having many cir- 
cumferential waves and has been referred to as plate-type flutter (ref. 20). The plate- 
type flutter is affected little theoretically by in-plane inertias and structural damping 
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(ref. 20). Thus a Donnell type theory, which should be sufficiently accurate for moderate 
to large values of 77 (ref. 22), is used and in-plane inertias and structural damping are 
neglected. 

Recent investigations have revealed that the aerodynamic boundary layer does not 
play an important role in shell flutter with many circumferential waves (ref. 21); thus, 
inviscid theory is used. Furthermore, it has been shown (ref. 21) that piston theory 
aerodynamics provides a convenient and useful tool for cylinder flutter and hence is used 
herein. 

A recent nonlinear analysis (ref. 23) indicates that cylindrical shell flutter can 
occur at aerodynamic pressures below the stability boundary for infinitesimal distur- 
bances. However, to make the present overall investigation tractable the small displace- 
ment infinitesimal linear flutter theory is employed. The stresses resulting from the 
loads p and P are obtained from the linear membrane solution, because stress dis- 
tributions from the more accurate solutions have been shown to have only a slight effect 
for a wide class of buckling problems. (See, for example, ref. 24.) Furthermore, the 
effects of displacements which occur prior to flutter, vibration, or buckling are neglected 
herein since they do not significantly influence cone buckling induced by lateral pressure 
(ref. 25) or the flutter characteristics of cylinders (ref. 26). 

Rings and/or stringers are assumed to be closely spaced so that their effect may 
be smeared out by use of equivalent orthotropic stiffness properties. Stiffness proper- 
ties were calculated from the formulas given in table I. The approach used herein 
neglects the effects of stiffener eccentricity. 

Governing Equations 

The principle of virtual work is used to derive the general nonlinear equilibrium 
equations and the linearized stability equations. A stress function is introduced so that 
the number of unknowns is reduced from three displacements to the stress function and 
the normal displacement. The introduction of the stress function necessitates the use of 
a compatibility equation. 

The equations are solved by an assumed displacement method. The normal dis- 
placement is assumed, the compatibility equation is solved exactly in terms of the arbi- 
trary coefficients in the assumed displacement, and then the equation for the displacement 
is solved by the generalized Galerkin method (refs. 15 and 16). Details of the derivation 
and solution of the governing equations are given in appendix B. These governing equa- 
tions are as follows: 


(y^yB^y + N y6B^ - n 0B = 0 

(2a) 

n 0B >( £ + (y N y0B) >y + N y0B = 0 

(2b) 
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D Y L i( w b) - N 0B c °t « - y N yA w B >yy - N 0 a| ^ + w B Jj - yp s = 0 (2c) 

where L^( ) is an operator defined by equation (B15). 

The boundary conditions on an edge y = Constant are to prescribe 


N yB - f ub = 0 

or 

U B 

(3a) 

N y6B " F vB = 0 

or 

V B 

(3b) 

QyB " F wB “ 0 

or 

W B 

(3c) 

MyB " M wB = 0 

or 

WB ,y 

(3d) 


Equations (2a) and (2b) are identically satisfied if the stress function F is introduced 
(refs. 27 and 28) so that, for rj > 0, 


n 0B _ F ,yy 

(4a) 

N yB = y F ,y + ^2 F ,00 

(4b) 

N y<® - F ,0 - f F ,y<P 

(4c) 

N 6>B = Fo^y 

(4d) 

NyB = jf 0 

(4e) 

N y0B = 0 

(4f) 


The stress functions are obtained from the compatibility equations (eqs. (B31) and (B34)). 
In terms of F, the compatibility equation is 

L 2 (F) = -Bg(l - MeMy)w B ^ yy cot o' (5a) 

where L 2 ( ) is an operator defined by equation (B33). For rj = 0 the equation for the 
stress function is 


F °,yy + y F °,y " By y 2° " "y M 1 " M e M y) w B >y cot “ 


(5b) 


The boundary conditions for the stress function are the boundary conditions for u B and 
v B , that is, equations (3a) and (3b). 
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To investigate flutter it is necessary that the surface loading p represent the 

s 

inertial and aerodynamic loading acting on the shell. If the piston theory aerodynamic 
approximation is used to represent the aerodynamic forces (for example, ref. 29), p s 
can be expressed as 

p s = -^tt - f *B >y - P c air w B )t <«> 

where in the denominator of the second term on the right the Mach number M^ has been 

replaced by /3 = \Jm^ - 1. Preference for the compressibility factor 0 is based on its 
demonstrated usefulness for the flutter of flat plates. 

With the use of equations (4a) and (6), the governing equation for the normal dis- 
placement wg (eq. (2c)), becomes 

D y L l( w B) - y N yA w B i yy - n 6A v B tU + w B >y ) - cot “ F ,yy 

+ y (yw B;tt + J- »B, y + P c alr w B ;t ) - 0 (’) 

For r] = 0, the term F V v is replaced by F 0 „. 

,y 

The flutter, vibration, and buckling characteristics of a truncated orthotropic coni- 
cal shell are now completely described by the equations for F (eqs. (5)) and the equa- 
tion for Wg (eq. (7)) with the appropriate boundary conditions (eqs. (3)). If the ortho- 
tropic character of the shell can be described by only distinguishing between Young's 
modulus and Poisson’s ratio in the y- and 0-directions, equations (5a) and (7) reduce to 
those obtained by Singer (ref. 30) for p s = 0. If the shell is isotropic, the equations 
reduce to those obtained by Seide (refs. 27 and 28) for p s = 0. 

The governing equations are solved in appendix B. The normal displacement Wg 
is assumed, and the equations for the stress function are solved exactly in terms of the 
arbitrary coefficients in the assumed displacement. The membrane stresses and 

NyA are related to the loadings p and P by the linear membrane expressions 
(eqs. (B36)), and the forces and moments applied to the shell edge by the elastic medium 
are expressed in terms of generalized stiffness and mass coefficients and shell displace- 
ments (eq. (B46)). The governing equation for wg is then solved by the generalized 
Galerkin method (refs. 15 and 16), which yields 
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1 1 



where a Q , a Q , and a m are the coefficients of the assumed normal displacement Wg, 
and 



The matrix elements by, Cy, dy , ey, and fy are given in the supplement to NASA 
TN D-5759. 

For buckling, the parameters C w and X are zero, either C p or Cp is speci- 
fied, and the remaining parameter becomes the eigenvalue. For free vibrations, Cp 
and Cp are specified, X is zero, and C w is the eigenvalue. For flutter, C w is 
again the eigenvalue but in addition to specifying C p and Cp, X and g a must also 
be specified. For g a > 0 it is necessary to specify the cone material, altitude (p and 
c air), and Mach number. For buckling and flutter, the circumferential wave number 17 
must be varied to determine the stability boundary. 
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For flutter, solution of equation (8) requires the extraction of complex eigenvalues 
for C w which occur as conjugate pairs. Thus, 


w = w Re + iw Im 

and 


( 10 ) 


C w ~ C cuRe + lC u>Im 


(ID 


The complex eigenvalues were found with the QR transform of Francis (ref. 31). As 
previously stated in the section "Assumptions," only the two lowest natural frequencies 
were examined to determine flutter. For g a = 0, flutter is defined as the lowest value 
of X for which C w becomes complex and the first two natural frequencies coalesce. 
For g a > 0, the aerodynamic damping coefficient g a is related to and C^ge 

by (see, for example, ref. 32) 


l C a>Iml 

tf c wRe 


(12) 


Further, g a can be related to X from the defining relationships given in equations (9); 
for isotropic shells the relation is 


^a 


0.662 


'm 


MT 


a /3 


(l - M 2 ) 


2 n 1/6 lv7h c air 


, 2/3 


M4 


(13) 


Thus, for g a > 0 flutter is defined as the value of X for which equations (12) and (13) 
give identical values of g a . 


RESULTS AND DISCUSSION 

Results are presented to indicate the flutter characteristics of simply supported 
isotropic conical shells. Sufficient numerical results are presented in the form of general 
flutter boundaries to permit the determination of the flutter condition for a wide range of 
cone angle, length-radius ratio, and radius -thickness ratio. In addition, results are pre- 
sented to indicate the trends of the effects of hydrostatic pressure, axial load, and ring 
and stringer stiffeners. Finally, results are presented to indicate the trends of the effects 
of finite elastic edge restraint provided by massless uncoupled springs; the limiting cases 
of free and clamped edges are also considered. Springs resisting only the meridional 
displacement Ug and the normal displacement wg (or the axial and radial displace- 
ments uy B and Ugg), the circumferential displacement v B , and the rotation w B 
are considered herein. In this investigation "conical" simple support is defined by 
Ny B = Vg = Wg = Myg = 0 and "cylindrical" simple support is defined by 
Vg = vg = Ugg = Myg = 0; clamped support is defined by ug = Vg = Wg = Wg = 0. 

9 J 
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Buckling and vibration results are presented in appendixes C and D. These results 
are compared with results of other investigations to indicate the accuracy and limitations 
of the present analysis. In addition, the trends of the effects of conicity, elastic edge 
restraint, applied loads, and ring and stringer stiffeners are indicated. Up to 40 terms 
were used in the computations to insure converged results. Results were assumed to be 
converged when the results (eigenvalues) for 1 + 4 terms differed by less than 1 percent 
from the results for I terms. 


Simply Supported Shells 

Effects of aerodyna mic damping.- The variation of the flutter parameter X with 
wave number ri for an isotropic conical shell with cylindrical simple support is shown 
in figure 2 along with the shell properties. The solid curve is for g a = 0; typical sta- 
bilizing effects of aerodynamic damping are indicated by the dashed curve. The minimum 
point of the flutter boundary was increased by about 8 percent for g a based on aluminum 
at sea level at Mi = 3. For denser materials or higher altitudes the effect would be even 
less. Calculations for cylinders based on the present theory (ref. 33) indicate the effects 
of aerodynamic damping on flutter increase with L/R and R/h. The magnitude of the 
effect was of the order of 20 percent for very large values of L/R and R/h but for 
practical values of these ratios is generally less than 10 percent. The effects of aero- 
dynamic damping are not considered in the remaining calculations. 

Effects of conicity .- Some effects of conicity on the flutter of isotropic conical shells 
with cylindrical simple support are shown in figure 3. Figure 3(a) shows the variation of 
X with cone angle a for two values of the ratio of length to radius at the small end of 
the cone L/Rj. The dashed portion of the curves represent extrapolation to a = 0°. 

The extrapolated result for l/r^ = 6 is in good agreement with the result of a corre- 
sponding cylinder presented in reference 34. For small values of a (a < 15°), X 
increases with increasing a. The phenomenon of increase in flutter resistance with 
increasing a has been obtained previously for three layered conical frustums (ref. 35). 
This trend is in marked contrast to the trends for buckling and vibration (appendixes C 
and D); the buckling load and natural frequencies always decreased with increasing a. 
However, the variation of X with L/Rj (fig. 3(b)) is similar to the trends for buckling 
and vibration. 

The increase in X with a is attributed to the fact that the shell flutter character- 
istics depend on the natural frequency spectrum as well as the magnitude of the frequen- 
cies. For example, a two-mode analysis of cylinder flutter (ref. 20) reveals that the 
critical value of the flutter parameter is proportional to A, the normalized difference in 
the square of the first two in vaccuo natural frequencies for a given value of 77 . Fig- 
ure 4 gives the variation of A with a for l/Rj = 6 for rj = 6 to 9. For a = 0° 
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the value of 77 for which A is a minimum (77 = 8) coincides with the value of 77 for 
which X is a minimum (table 5.2 of ref. 34). For a > 0° the values of 77 for which 
A and A were minimums differed only slightly; the difference tended to increase with 
increasing a . Thus, A may be useful for predicting the critical value of wave number 
for flutter of both cylindrical and conical shells. For small values of a (a < 15°), the 
effect of increases in frequency spacing with a is apparently greater than the effect 
of decreases in frequency magnitude and A increases as does the flutter resistance. 
However, for large values of a (a > 15°), increases in a affect the frequency spacing 
and magnitudes in such a manner that A decreases, as does the flutter resistance. 

The variation with L/Rj of the ratio of the maximum value of q//3 to the value 

for a ss 0° is indicated in figure 5. As can be seen the variation of — with 

L /R^ is essentially linear and depends on the value of Ri/h. For R^/h = 200, 

(q//3)max occurs at a= 0° for L/R^ S 1. 

For a cylinder, A is based on free-stream conditions, which are the same as local 
conditions, whereas for a cone, A is based on local conditions behind the cone bow shock. 
In order to determine the trends for A based on free-stream conditions, the results for 
L/Rj = 6 were converted to free-stream conditions by use of reference 36. The free- 
stream Mach number Moo was taken to be 3.0. The results for both local and free- 
stream conditions are shown in figure 6 in terms of A and a. Comparison of these 
results reveals that use of free-stream conditions can significantly reduce the magnitude 
of A for a; > 5° but does not alter the overall trends. 

Comparison with previous in vestigation .- Shulman (ref. 34) presents flutter results 
for a simply supported conical shell with L/Rj = 8.13, Ri/h = 148, and a = 5°. These 
proportions are such that 1 / R av = 6 and R aV/ /h = 200. In the reference, results for 
the cone were compared with results for a cylinder with L/R = 6 and R/h = 200 to 
indicate the effects of conicity on flutter. The flutter analysis was formulated on the 
basis of Donnell shell theory and piston theory aerodynamics. The stress function was 
eliminated from the governing equation for the displacement w, and the resulting eighth- 
order equation was solved by the Galerkin method. It should be noted that not all bound- 
ary conditions were satisfied by the assumed displacement for the cone. 

The results of reference 34 for eight and four terms are shown in the following 
tables in terms of the parameter of the present investigation along with the results from 
the present investigation for the same two shells; for both investigations, g a = 0. The 
results for the cylinder from the present investigation were obtained by extrapolation. 

For L/R av = 6 and R av /h = 200 and for a = 0°, 
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A for — 


V 

Present investigation 

Reference 34 


8 terms 

8 terms 

6 

605 

605 

7 

410 

412 

8 

355 

352 

9 

370 

360 

10 

425 

427 


and for a = 5° 





X for - 

- 


V 


Present investigation 


Reference 34 


4 terms 

6 terms 

8 terms 

12 terms 

4 terms 

4 

812 

1037 

1056 

1074 

1251 

5 

492 

577 

588 

590 

688 

6 

622 

582 

604 

607 

669 

7 

1042 

] 

627 

643 

652 

714 


Examination of these tables reveals that the results of the present investigation are 
in excellent agreement with the results of reference 34 for the cylinder. The results of 
reference 34 for the 5° cone indicate that conicity has a large stabilizing effect, the criti- 
cal value of A being increased by nearly a factor of 2. The results of the present inves- 
tigation also indicate a similar stabilizing effect of conicity but slightly less than the 
increase given by the results of reference 34. In the flutter parameter used in refer- 
ence 34, there appears a quantity defined by equation (3-31e) of reference 34 to be 

f(6) = In 6* = i tanh -1 6 

where 


6 * Ri 


r 2 1 + 6 


1-6 


There appears to be a typographical error in the first definition and it should read 

f ^=^=| tanh ' 16 

The reasons for the larger difference in flutter results for the 5° cone as compared 
with the flutter results for the cylinder are attributed in part to the use in reference 34 of 
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the Galerkin method in conjunction with the uncoupled eighth-order shell equation and an 
assumed displacement that did not satisfy all boundary conditions; it has been demon- 
strated for buckling of cylinders (refs. 37 and 38) that application of the Galerkin method 
to the governing eighth-order equation on the normal displacement can give incorrect 
results for the buckling load or mode. Some comments on why application of the Galerkin 
method to the uncoupled equation can lead to erroneous results have been given by Singer 
(ref. 39). 


Flut ter mode shapes .- The effects of cone angle on flutter mode shape is shown in 
figure 7 for various values of a for isotropic conical shells with cylindrical simple 
support with L/R^ = 2 and Rj/h = 200. The maximum amplitude occurs near the 
trailing edge, and the position of maximum amplitude moves towards the trailing edge as 
a increases. The mode shapes are similar to those calculated for flat plates (ref. 29) 
but differ from the one calculated for a cylinder with L/R = 1.93 (ref. 21) in that the 


cylinder flutter mode did not have a node line away from the boundaries. Figure 8 shows 
the effect of length-radius ratio L /R\ on flutter mode shape of isotropic conical shells 


with cylindrical simple support with Ri/h = 200 and a = 30°. For large values of 
L/Ri (l /Ri = 6j, several node lines occur in the flutter mode as indicated in figure 8(c). 


Genera l flutter boundaries .- The variation of X with Rj/h for conical shells 
with cylindrical simple supports was found to be well correlated by the parameter X 


determined in reference 33 for cylinders, where, for conical shells, 


X = — — where x = tanh 

(i-m 2 )(V>) x 

The variation of X with L/R^ is shown in figure 9 for several values of a. Super- 
sonic flow over a cone in air at zero angle of attack is dependent on cone angle a and 
free -stream Mach number M m . For example, at = 3, a must be less than 
about 46° for > 1 (ref. 36), and for = °°, a must be less than about 57.5°. 
However, supersonic flow can occur over part of a cone for a > 57.5° if the cone is at 
an angle of attack. Since 60° cones are of interest as aerodynamic decelerators (ref. 25) 
which can be at angles of attack, results are presented in figure 9 for 0° lal 60°; the 
results for a = 0° were obtained from reference 33. The curves represent average 
values of X obtained for values of R,/h from 100 to 2000. No calculated value of X 
differed by more than about 12 percent from the average values shown in figure 9. Thus, 
these curves provide a rapid and reasonably accurate means of predicting the flutter char- 
acteristics of isotropic conical shells with cylindrical simple support for a wide range 
of a, L/Rj, and Rj/h. Blunt cones are significantly more susceptible to flutter than 
cylindrical shells for the same local flow conditions when the radius of the cylinder is 
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equal to R^. Results for L/R^ = 1 could not be correlated by use of the parameter A. 
Difficulty in correlating results in this region might have been expected from the trends 

(q /^max / , 

exhibited by the variation of — — — — with Rj/h for L/Rj S l (fig. 5). 


(q/ 0)a=O 


Effects of ap plied load.- Some effects of hydrostatic pressure and axial load on the 
flutter of isotropic conical shells with cylindrical simple support are shown in figure 10 
in terms of A as a function of the ratios p/p cr (positive for external pressure) or 
P/P c r (positive for compression) for a = 10° and 30°. The variation of A with applied 
load is nearly linear for both hydrostatic pressure and axial load. The effect of hydro- 
static pressure is slight, whereas the effect of axial load is pronounced. Compressive 
axial load can decrease A by up to an order of magnitude from the value for no load; 
similar trends have been shown for cylinders (refs. 40 and 41). Thus the effect of Ny^ 
(the only stress for axial load) is considerably greater than the effect of N (the pre- 
dominate stress for hydrostatic pressure). Results for flat plates also indicate that the 
effect of the stress in the direction of the airflow is much more significant than the stress 
in the cross-flow direction (ref. 42). 


Effects of compressive axial load on flutter mode shape are shown in figure 11 for a 
conical shell with cylindrical simple support for a = 10°. As can be seen, the effect of 
axial load is to move the point of maximum amplitude toward the leading edge and to 
increase the amplitude over the forward part of the shell. Similar effects of compressive 
load have been observed both theoretically and experimentally for flat plates (ref. 43). 


Effects of r ing or stringer stiffening.- In order to analyze a ring- or stringer- 
stiffened cone by use of the present analysis, it is necessary to calculate equivalent ortho- 
tropic stiffness properties. Stiffness properties which neglect eccentricity effects were 
calculated from the formulas given in table I. Some flutter trends for ring- and stringer- 
stiffened shells are shown in figure 12 for a = 30° with stiffener details also given. 
Stiffener size and spacing was such that the ring- or stringer-stiffened shells were essen- 
tially of equal weight. The stringer spacing, which actually varies along the length of the 
shell, was assumed to be the constant value of the spacing at the small end of the shell in 
calculating the equivalent orthotropic stiffness constants. Buckling and vibration results 
for the same shells are presented in appendixes C and D for 5° i a S 60° and in refer- 
ence 44 for a = 0. 


The flutter results are presented in figure 12 in terms of the ratio of the value of 
q//3 for the stiffened cone to the value of q/j3 for the unstiffened cone plotted against 
L/Rj. The addition of either rings or stringers for the present configuration increases 
the value of q//3 required for flutter compared to the value for the unstiffened shell. 

The effect of stringer stiffening decreases rapidly with L/Rj up to about 3 and becomes 
essentially independent of l/R^ for L/Rj > 5. On the other hand, the effect of ring 
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stiffening increases with L for the entire range of L/Rj. The results presented 
in figure 12 suggest for the present configuration that, for L/Rj less than about 1, 
stringers may have more effect than rings. The effect of both rings and stringers 
decreases with increasing cone angle as shown in figure 13; the reduction is more pro- 
nounced for rings. The addition of stringers has a larger effect on flutter than on either 
buckling or vibration. This result is attributed to the fact that the effect of stringers 
increases with the number of axial half-waves m, and vibration modes with m > 1 have 
a pronounced effect on flutter characteristics. 

The results presented in figures 12 and 13 neglect the effects of stiffener eccentric- 
ity. Two mode results for the flutter of flat plates with the same stiffener details shown 
in figure 12 are presented in reference 44. These results indicated that consideration of 
stiffener eccentricity increased the flutter resistance by nearly a factor of 2. Thus it 
would appear that stiffener eccentricity effects would be significant for the flutter of shells 
also, although there are no results available to verify this conjecture. 


Elastically Supported Isotropic Shells 

Effects of axial and rotational restraint .- Effects of axial restraint (kj) are shown 
in figure 14, which gives the variation of X and A with the ratio kj^E; shell proper- 
ties are also given. Increasing axial restraint results in decreases in X of up to about 
10 percent for kj/E = 1. This result is in marked contrast to the trends for buckling 
(appendix C) and vibration (appendix D); the buckling load and natural frequencies always 
increased with increasing kj/E. The decrease in X with kj/E is attributed to the 
dependence of shell flutter characteristics on the frequency spectrum as well as magni- 
tude. The decrease of A with k^E indicates that the effect of decreases in frequency 
spacing with k^/E is greater than the effect of increases in frequency magnitude. These 
results illustrate again that trends established for buckling and vibration do not necessar- 
ily hold for flutter. 


Some effects of rotational restraint are shown in figure 15, in terms of X and 

Rik 4 

— - — for the same shell considered in figure 14. The lower curve is for complete axial 

D y 

constraint, the upper curve for no axial restraint. Figure 15 shows that X increases 
only slightly with increasing rotational restraint. Calculations (not presented) revealed 
that A also increased with increasing rotational restraint. The increase of X with 
R l k 4 

— is similar to the trends for buckling and vibration. The results also indicate that 

D y 

the value of X for the essentially clamped cone is about 6 percent less than the value of 
X for the simply supported cone. However, additional calculations (not presented) indi- 
cated that the effects of edge clamping depend on cone angle a. As a approached zero 
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(a < 15°), the effects of edge clamping (primarily axial restraint) became stabilizing as 
for cylindrical shells (ref. 26). 

Effects of circumferential and radial restraint .- Effects of elastic circumferential 
restraint k£ and radial restraint kg are shown in figure 16, which gives the variation 
of X with the ratio kg/E; shell properties are also given. The upper curve is for 
kg/E = 1 (vg = 0), the lower curve for kg/E = kg/E; for both curves kj = k 4 = 0. The 
numbers on the curves represent the number of circumferential waves in the flutter mode. 
The open circles and squares at kg/E = 1 (ugg = o) show the slight effect of changing 
the circumferential restraint from complete constraint (vg = o) to no restraint (Nygg = o) 
for no radial displacement. The effect of decreasing the radial restraint becomes signifi- 
cant for kg/E i 10"3. For kg/E = 1, decreasing kg/E to 0 decreases X by over a 
factor of 2 as indicated by the solid square. However, when kg/E = kg/E, decreasing 
kg/E to 0 decreases X by over 2 orders of magnitude as indicated by the solid circle, 
and the flutter mode changes from 77 = 8 to an essentially inextensional mode at 77 = 2. 
The reason for the large effect of circumferential restraint as the radial restraint 
becomes weak is that the constraint Vg = 0 prevents the inextensional form of flutter 
(77 = 2), which is the critical mode of instability when all edge restraints are relatively 
weak. The flutter trends shown in figure 16 are similar to the trends for free vibration 
for the same shell given in appendix D. Since the critical value of 77 becomes small as 
the edge restraint becomes relatively weak, the present Donnell type theory becomes more 
approximate. Hence it might be desirable to use more exact shell theories in analysis of 
flutter of shells with weak edge restraint. 

Calculations for the variation of A with the various edge restraints considered in 
figures 14, 15, and 16 indicated that the qualitative effects of edge restraint on X appar- 
ently can be determined by merely examining the effect of edge restraint on A for the 
appropriate range of 77 . 


CONCLUDING REMARKS 

The flutter, vibration, and buckling characteristics of orthotropic truncated conical 
shells with generalized elastic edge restraint have been investigated theoretically. The 
problem was simplified by making certain assumptions which have been found admissible 
for flutter analyses of cylinders. The shell analysis utilized the classical Donnell type 
theory, in-plane inertias and structural damping were neglected, and the aerodynamic 
loading was represented by the inviscid two-dimensional quasi-steady approximation 
(modified piston theory). The principle of virtual work was utilized to formulate the 
problem. An approximate solution was obtained by the generalized Galerkin method in 
terms of assumed displacement functions that had to satisfy only geometric constraints. 
A stress function was introduced so that only the normal displacement needed to be 
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assumed. The equation for the stress function was solved exactly in terms of the arbi- 
trary coefficients in the assumed displacement. 

Numerical results for flutter of conical shells revealed that trends established for 
buckling and vibration do not necessarily hold for flutter. Increasing cone semivertex 
angle (up to about 15°) for length-radius ratios greater than about 1 increased the resis- 
tance to flutter even though the buckling load and natural frequencies decreased. Also, 
increasing the axial edge restraint could decrease the resistance to flutter even though 
both the buckling load and natural frequencies were increased. These phenomena were 
attributed to the fact that flutter depends on frequency spectrum as well as magnitude. 

The flutter resistance for a given circumferential wave number was shown to be propor- 
tional to the difference in the square of the first two in vaccuo natural frequencies. Thus 
the effect of a given parameter on the flutter of a conical shell can be determined quali- 
tatively by considering its effect on the natural frequencies (magnitude and spectrum) of 
the shell for the appropriate range of wave numbers. Although slightly conical shells 
(semivertex angles less than about 15°) have more flutter resistance than cylinders with 
the same thickness and length and a radius equal to the smaller radius of the conical shell, 
blunt cones (semivertex angles greater than about 45°) are significantly more susceptible 
to flutter than cylinders. Sufficient numerical results are presented in the form of gen- 
eral flutter boundaries to permit the determination of the flutter condition for simply 
supported isotropic conical shells for a wide range of cone semivertex angle, length- 
radius ratio, and radius -thickness ratio. Results are also presented to indicate some 
effects of finite axial, radial, circumferential, and rotational edge restraints, hydrostatic 
pressure, axial load, and rings and stringers on the flutter, vibration, and buckling char- 
acteristics of conical shells. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 18, 1970. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Con- 
ference on Weights and Measures, Paris, October 1960 (ref. 45). Conversion factors for 
the units used herein are given in the following table: 





— -- 

Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

SI Unit 



(*) 


Area 

in 2 

0.6451 X 10“ 3 

meters^ (m^) 

Length 

in. 

0.0254 

meters (m) 

Moment of inertia 

in* 

0.4161 X 10~ 6 

meters^ (m^) 

Young's and shear moduli . . . 

lb/ in^ 

6.895 X 10 3 

newtons/meter 2 (N/m^) 


* 

Multiply value given in U.S. Customary Unit by conversion factor to obtain 
equivalent value in SI unit. 


Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

giga (G) 

10 9 

kilo (k) 

10 3 

centi (c) 

1— 1 
O 
l 

CO 

milli (m) 

10" 3 
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APPENDIX B 


GOVERNING EQUATIONS 

The principle of virtual work is used to formulate the problem of flutter of an 
orthotropic conical shell with generalized elastic edge restraint. Such a formulation 
admits a solution by the generalized Galerkin method in terms of assumed displacements 
that must satisfy only geometric constraints. For a generalized elastic edge restraint 
there are no geometric constraints (only force or natural constraints); thus a rather 
simple series can be used for the assumed displacements. Such a procedure is straight- 
forward and well known (for example, refs. 15 and 16) but has been used only sparingly 
(for example, refs. 46 and 47). Instead, investigators of flutter problems have obtained 
displacement functions that satisfy all boundary conditions (sometimes going to great 
lengths, for example, ref. 48), or have relied on psuedo energy methods (for example, 
refs. 32 and 49). 

Although the use of complicated series representations can result in fewer terms 
required for convergence, and thus possibly justify the effort, the modification of potential 
energy methods for nonconservative problems seems unjustified. As pointed out in ref- 
erences 39 and 50, investigators have applied the Galerkin method to single uncoupled 
equations for systems properly defined by two or more coupled equations and obtained 
erroneous results. The erroneous results are due to inadequate attention to the boundary 
conditions that result from the manipulations required to obtain the single uncoupled equa- 
tion. Thus in the present investigation the derivation and solution of equations is given 
in some detail to emphasize the straightforward yet completely general character of the 
virtual work (generalized Galerkin method) approach to nonconservative problems. 


General Nonlinear Equilibrium Equations 


Donnell type expressions for middle surface strains and curvatures of a deformed 
conical shell are given by Seide in reference 27. For this paper, these expressions are 
given in the following form: 


e y = u ,y + ?( w ,yf 




€ e " 


u - w cot a 1 1 / \ 

y y ,<P iv2( »^) 


2y <! 


6 ye = V ‘ y + ^ + f w .y w ,0 




(Bl) 
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*y = w ,yy 

Xe = | w ,y + ^2 w ,4>4> > 
x y0 _ y w >y<t> ~ Z2 w >0 

Middle-surface force and moment resultants are defined as 

N y = B y( e y + M e 6 e) 

Ng — Bg^eg + M y^yj ^ 

Nyg = BygCyg 


(B2) 


(B3) 


My = -Dy( Xy + MgXg) 
Mg = -Dg^Xg + MyXy) /> 


M y 9 ~ - D y0X y g 


(B4) 


The equations of equilibrium will be obtained by the principle of virtual work, which 
requires that 

Sfi t = SU s (B5) 


where 60* represents the total virtual work of all the external forces and the restoring 
forces of the elastic edge restraint, and 6U S represents the change in the internal strain 
energy of the system corresponding to the virtual displacements. 

The virtual work of the total surface pressure can be expressed by the surface 
integral 

n y 9 p27rsino! 

6 Op = J J p t 6w y dcp dy (B6) 

yi ^ 


while the virtual work of the boundary forces can be expressed as 
>27TSinQ!/ 

'0 


p27TSinQ!/_ _ _ 6 w A 

60j = j (N y 6u + N y g6v + Q y 6w - My6w^y + Myg — — |y 


y2 


d<£ 


(B7) 


yi 


where 6 represents virtual change in variable and (p = 8 sin a. The effects of the 
generalized elastic edge restraint can be introduced into the analysis by the virtual work 
expression 
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n 2jr sin a . . 

= ^F u 6u + F v 6v + F w 6w “ M w 8w ? yjy 


y2 

yi 


d 0 


(B8) 


where F U j F V> F Wj and M w are the forces and bending moment exerted on the shell 
by the elastic support at the shell edge. The total virtual work of all the external forces 
is given by 

8f2^. = 8&2p "i* 80^ + 8£2g (B9) 


The strain energy of the shell can be expressed in terms of middle surface force 
and moment resultants by the surface integral 


U s = 



+ N 0 e e +N ye e y g - M yXy - M 0X0 - 2M y0 X y0 )y d 0 d Y 


(BIO) 


If the force and moment resultants, strains, and curvatures are expressed in terms of 
displacements by equations (Bl) to (B4), if the variation of Us is taken with respect to 
these displacements, and if equation (B9) is utilized, equation (B5) becomes, after integra- 
tion by parts, 



- jeq. (B16d)j Sw >y j y 


y2 


d0 = 0 


yi 


(Bll) 


where the general nonlinear equilibrium equations are given by 


( yN y),y + Ny0 ,0 ■ N<? = 0 

(B12) 

Ne ,0 + ( yN ye),y +N ye = ° 

(B13) 
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(B14) 
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with 


Lj( ) 


: vt ) 


^ ^ \y<P<P + y^ \yy<P<P 

and the boundary conditions on an edge y = Constant are to prescribe 


% - F U - Ny = 0 

or 

U 

Ny0 - Fy - Nye = 0 

or 

V 

Qy - F w - Qy + ^Mye)^ = 0 

or 

W 

My - M W - My = 0 

or 

w 

N fl w . M fl ( 
% = N y w ,y + - -jr + - 

yM y),y + 
y 

2M y«, 0 

y 


+ 2 


Dy + " 6 


(B15) 


(B16a) 

(B16b) 

(B16c) 

(B16d) 

(B17) 


Linearized Stability Equations 

The total surface pressure loading p^. can be expressed as 

P t = P + P s (B18) 

where p represents an axisymmetric uniform lateral pressure loading that exists prior 
to flutter and p g represents the asymmetric surface loading resulting from vibration 
and flutter of the shell. Because the loading prior to flutter is axisymmetric 


NyeA = V A = ( = 0 


(B19) 


where the subscript A denotes quantities that exist prior to flutter, vibration, or 
buckling. If equations (B19) are utilized, equation (B13) is identically satisfied and equa- 
tions (B12) and (B14) become 


D, 


yw» + 2w 

A ,yyyy 


(yN yA ),y - n„a = 0 
D 


(B20a) 


A yyy + D 


ei 


y\y 


,2 W A 


*y 


-? W A 


yy 


- Ng A cot a 


- (yN yA w A;y ) y - yp - o 


(B20b) 
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and the boundary conditions on an edge y = Constant are to prescribe 


S5 

i 

> 

i 

55 1 
<< 

il 

o 

or 

U A 

(B21a) 

> 

ti 

o 



(B21b) 

QyA ■ F wA - Qy = 0 

or 

w A 

(B21c) 

M yA " M wA - M y = 0 

or 

W A,y 

(B21d) 


If the exact nonlinear stress distribution prior to flutter were desired, it would be neces- 
sary to solve equations (B20) subject to the boundary conditions given by equations (B21). 


To investigate flutter, the total displacements are expressed as the sum of the axi- 
symmetric displacements that exist prior to flutter and the infinitesimal asymmetric dis- 
placements ug, Vg, and Wg that occur at flutter. Thus 


u = u a + u b 


"\ 


V = V T 


w = w A + w B 


(B22) 


Similar expressions can be developed for the stress resultants, moments, and forces 
applied by the elastic edge support. 

If these expressions and equations (B18) and (B22) are substituted into equation (Bll), 
if equations (B19), (B20), and (B21) relating subscript A terms are subtracted out, if 
products of the displacement w^ and the stresses Nyg and Nygg are omitted, and 
if terms nonlinear with respect to subscript B are neglected, the virtual work expres- 
sion governing stability of the conical shell becomes 



(B23) 
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where the linear stability equations are 

(yN yB ) jy + N y 6> Bj 0 - N eB = 0 (B24) 

N 0B,0. + ( yN y0B) >y + N y 6 B = 0 ( B25 ) 

D y L i( w B) " N eB cot 01 ~ yN yA w B ?yy " n 0A^f + w B y ) - yp s = 0 ( B26 ) 


and the boundary conditions on an edge 

N yB " F uB = 0 
N y0B " F vB = 0 
QyB “ F wB = 0 

M yB " M wB = 0 
In obtaining equation (B26), use has been 


Constant 

are to prescribe 


or 

U B 

(B27a) 

or 

V B 

(B27b) 

or 

W B 

(B27c) 

or 

WB ,y 

(B27d) 


made of the identity 


( yN yA w Bj y) y + |( N eAWB i0 ) >0 = y N yA w B t yy + Ne(i Wg ^ + w B , 

+ W B,y ( yN y A ),y “ N 0A 


(B28) 


From equation (B20a), the term in brackets in equation (B28) is 0. 

Equations (B24) and (B25) are identically satisfied if the stress functions F or 
F 0 (refs. 27 and 28) are introduced such that, for 77 > 0, 


n 0B = F ,yy 

N yB = 7 F y + ^ F 00 > 

NyeB = F >y < t> 


(B29) 
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and, for 77 = 0 , 


N 0B = F c 


F 

N y B = ir 


Ny0 B = 0 


-\ 

,y 




j 


(B30) 


where 77 is the number of circumferential waves. 

Thus, the number of unknowns is reduced from three (ug, Vg, and Wgj to two 
(F and Wg). However, in order to use the stress -function approach an equation is 
needed to define F. The appropriate equation is the compatibility equation which is 
obtained by eliminating ug and Vg from the linear form of the middle-surface strain 
displacement relations and is given as follows: 

( y£ yeB,4 y - £ yB, w - (y%B, 4 * y*yB, y - T «* » "B, yy 

or, in terms of F, 

L 2 (F) = -B fl (l - d’ e My)w Bjyy cot a 


(B31) 

(B32) 


where 


Bz 


l 2 ( ) = y( ),yyyy + gT 


1 2 11 

),0000 + -A ),00 + 3< ),y “ y< >, 


yy 


^(l “ M>y) - V 0 ~( ) >( p (j) - ^ ),y00 ),yy <p<p 


(B33) 


The boundary conditions for F are the boundary conditions for ug and Vg 

(eqs. (B27a) and (B27b)). For 77 = 0, the symmetry conditions N y0B = Vg = ( = 0 

exist and the compatibility equation becomes 



-Wtj cot a 

B ,y 


(B34) 


or, in terms of F 0 , 


F °,yy + ? F °, 


l!S = -5 Be ( 1 -' J >y) cot “ WB 


By y 2 


»y 


(B35) 


The boundary condition for F 0 is the boundary condition for Ug (eq. (B27a)). The 
linear membrane stress distribution is adequate to express and N y A in terms of 
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the external hydrostatic pressure p and externally applied axial load P (positive in 
compression) . Thus , 


N eA = -y tan a p 


N vA = tan a p - : 

y A 2 27 t sin a cos o' 


(B36) 


To investigate flutter it is necessary that p s represent the inertial and aerody- 
namic loading. If the modified piston theory aerodynamic approximation is used to repre- 
sent the aerodynamic forces (ref. 29), p s can be expressed as 


p s = " rW B,tt 


2q 

J8 


w 


B, y - Pc air w B ;t 


(B37) 


where the first term on the right is the inertia loading, the second term is the static aero- 
dynamic loading, and the last term is a viscous loading corresponding to the aerodynamic 
damping. In the denominator of the second term on the right, the Mach number has 


been replaced by /3 


4 


m i - »■ 


Preference for the compressibility factor f3 is based 


on its demonstrated usefulness for flutter of flat panels. 


With the introduction of the stress function and equations (B36) and (B37), the virtual 
work stability criterion becomes 


pi 

J yi J o 


277 sin a 


n 2 tt sin a r 

Jeq. (B39)j 6w R d <p dy + J jeq. (B27c)J 


5w 


B 


- (eq. (B27d)J 6w B y 


y2 


dcf) = 0 


(B38) 


yi 


where the governing equation for the normal displacement w B is given by 

D y L l(w B ) - cot « F_ yy + (y ^ P + 2, - ain « cos a ) y"B,yy + * “ p(| + W B y ) 


+ y(rw B;tt + j- w B>y + PCair .w Bjt ) - 0 

For 7] = 0 the term F yy is replaced by F 0 y 


(B39) 
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Solution of Equations 

The equations for the stress functions (eqs. (B32) and (B35)) and the equation for 
the normal displacement (eq. (B38)) are solved by an assumed displacement method. 
Equation (B32) has variable coefficients that reduce to constant coefficients by use of the 
transformation of coordinates (ref. 51) 

y = y^e x (B40) 

y 2 

Thus, y = y^ corresponds to x = 0 and y = y^ corresponds to x = x^ = In — . For 
simple harmonic motion the displacement Wg is assumed to be 

Wg = Wg(x) cos ncp e iwt (B41) 


where 


M 

Wg(x) = a Q + a 0 x + ^ a m sin £2 
m=l 1 


(B42) 


and 


n = 


U 

sin a 


(B43) 


The expression for Wg is continuous in the circumferential direction and permits deflec- 
tions and rotations to exist at the ends of the cone. 


The boundary conditions (eqs. (B27)) are now specialized, 
displacements are of the form (neglecting the term e iCt> ^) 


The shell stresses and 


U B = u b cos n 0 
v B = Vg sin ncf) 

w B = Wg cos n 4> 

N yB = N yB COS n< ^ 

N 0B = N 0B COS n(f) 
N y6B = N y0B sin n t 


(B44) 
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From equations (B27) it follows that (the expressions for Q yB and M yB being 
considered) 


-\ 


F uB = F uB COS n(p 


F vB = F vB Sin n(p 


F wB = F wB C0S ncf> 




M 


wB = M wB COS ^ 


(B45) 


The forces F* B , F* B , F^ B , and M^ B are now expressed in terms of the shell dis- 
placements u B , v^, w B , and w^ . The most general relations for a shell of revolu- 

,y 

tion supported by some generalized unloaded elastic support at its boundaries can be 
written in matrix form as (ref. 17) 


uB 


vB 


wB 
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u^ 
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(B46) 


w T 




'Z, 


The term Kij represents the elastic restraint of the support at the shell edge, and the 
term represents the mass effects of the elastic support. By utilizing equa- 

tions (B44) to (B46), the boundary conditions (eqs. (B27)) become 


N 
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-M yB 
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U, 
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>y 




where the minus sign applies at y = y^. 
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(B48) 


(B49) 


where 

kjg = - kg^ cos a sin a (B50) 

and k^ and k^ are spring constants. 

The nondimensional forms of Kjj and G^j are Sy and My; expressions for 
Sjj and My relating the stiffness and inertia characteristics of end rings to the ring 
and shell material properties and geometry were obtained from equations (A5) to (A8) of 
reference 17 and are presented in the supplement to NASA TN D-5759. For flutter (or 
vibration), the eigenvalue u> appears in the boundary conditions unless all terms 

are 0; thus, it is necessary to iterate on the eigenvalue until the value used in the bound- 
ary conditions differs by a negligible amount from the calculated eigenvalue. 

C ompatibility equations .- The procedure for solving the compatibility equation for 
the stress functions is outlined in this section. The details of the solutions are given in 
the supplement to NASA TN D-5759. If the stress function F is taken to be 

F = e x f(x) cos n cp e ia * (B51) 

and equations (B40) and (B41) are utilized, equation (B32) reduces to the ordinary differen- 
tial equation 

f,xxxx + Pif jXX + P 2 f = -Bei 1 ' ^y) y l cos a ( W B ^ " W B x ) (B52) 
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where 




p 2 = |V-i) : 
B y 


b. 


-T\ 


B. 


ye 


(l - M^My) - 2^0 




(B53) 


By use of equation (B51), equation (B52) contains only even derivatives of f with 
respect to x. The homogeneous solution is readily found to be 


where 


AiX -Aix AoX -Aox 

f^ = A^e 1 + Age 1 + A^e + A^e 


(B54) 


X 1 = 


X 2 “ 


/ 2 \l /2 

-Pi + (p* - 4P 2 ) 

1/2 1 

2 

\ 

r 1 2 n 1 / 2 ' 

- p i - ( p i - 4p 2 ) 

1/2 

2 
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(B55) 


and Aj to A 4 are constants of integration. The particular solution for an M-term 
approximation in Wg is obtained by substituting 


fp - ^o 


M . 

It 


d m cos 


mnx 


+ e 


m 


sin 


mwx 


(B56) 


into equation (B52) and matching coefficients. With f = f^ + fp known in terms of Aj 
to A 4 , the stress resultants are obtained from equations (B29) and (B30), the strains are 
obtained from the constitutive equations (B3), and the displacements ug and vg are 
obtained from the linearized strain-displacement relations (eqs. (Bl)). The constants of 
integration A 4 to A 4 are determined from the boundary conditions for ug and vg, 
the first two of equations (B47). Expressions for the stress function, stress resultants, 
displacements, and constants of integration are given in the supplement to NASA 
TN D-5759. 

Stability equations.- By utilizing the transformation of coordinates (eq. (B40)) the 
virtual work expression for the stability of the conical shell becomes 
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r x 1 r> 2 it sin ex r 
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where the equation for Wg is given by 
cot a ( f ,xx + f ,x) 
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277^yje x ^ sin a cos a 
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(B58) 
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(B59) 


For 77 = 0 the term 


cot a (f >xx + f jX ) 


y?e x 


in equation (B58) is replaced by 


cot a F, 


y^ 


If F wB and IV^g in equations (B27c) and (B27d) are expressed in terms of Kjj, 
Gjj, and the shell displacements by use of equation (B46), if the moments, transverse 
shear, f, Ug, Vg, and Wg are expressed in terms of the assumed series for Wg, and 
if the integration with respect to <p is performed, equation (B57) becomes the general- 
ized Galerkin solution of equation (B58) where 5wg and 6wg are taken successively 
as 
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~\ 


5w^ = 6a 0 


6wg = Sa Q x 


6w B = 6a P sin 


5w* = 0 

B,x 


6w* = 6a 0 

B,x 




6Wg x = 5a p EE cog Hi (p = 1,2,3,.. .,M) 


1 X 1 


J 


If the integration with respect to x is also performed, equation (B57) becomes M 
simultaneous algebraic equations which can be written as 
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and b^, c^, 


djj, e^, and f^ are given in the supplement. 
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(B60) 


+ 2 
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BUCKLING RESULTS 
Simply Supported Shells 

Effects of conicity .- Some effects of conicity on the buckling of isotropic truncated 
conical shells with cylindrical simple support subjected to hydrostatic pressure are 
shown in figure 17. Figure 17(a) shows the variation of Cp with cone angle a for 
two values of L/Rj. The dashed portion of the curves represent extrapolation to zero 
cone angle. These results are in excellent agreement with the results obtained in ref- 
erence 44 indicated by the circles. As would be expected, the buckling load for a fixed 
small-end radius R 1 decreases considerably with increasing a; for a = 45° the 
buckling pressure is about an order of magnitude less than the pressure required to 
buckle a cylinder with the same L/Ri and R^/h. Figure 17(b) shows the decrease of 
Cp with increasing L/Rj for a = 0° and 30°. Again, the extrapolated results for 
a = 0 are in excellent agreement with the results of reference 44. 

Buckling results for isotropic conical shells have been correlated (ref. 52) by use 
Ri 

of the taper ratio 1 - =— where 

r 2 

Ro T . 

— = 1 + — tan a (Cl) 

Rl Ri ' ' 

and the concept of an equivalent cylinder. Calculated results for 3° § a £ 75°, 

T Ri 

1 S — — = 6.28, and 250 ^ — = 1000 are presented in figure 18 in terms of the ratio 
Ri h 

P R 1 

as a function of the taper ratio 1 - =— where p , is the buckling pressure for 

P C yl Rg Cyl 

an equivalent cylinder. If the normal equivalent cylinder transformations 


(C 2) 
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are substituted into the cylinder buckling equation (eq. (17) of ref. 53), the following 
expression for p cyl for unstiffened or ring-stiffened shells is obtained for buckling 
into one-half wave (m = 1): 


cyl 
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2? j_ sin^ a cos a E 
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( R 2/ R l) 
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(C3) 


(C4) 


It should be noted that equation (17) of reference 53 contains expressions multiplied 
by z R , which can be either positive or negative. These terms, which account for eccen- 
tricity (one-sided) effects, are omitted in equation (C3). 

Equation (C3) is a function of 77, which must be varied to find the minimum values 
of p C yp An approximate expression for p cy j, based in part on the assumption that 


[( R 2/ R l)+l] „ 24 

[( R 2 / R l) - ’ Si " “ 

is (ref. 54) 
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where for ring-stiffened shells 


D fl 

— - = 1+12 
D y 




e r 

J R 
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v 3 
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a R / z sR ^ 2 
d R h\ h / 


(C7) 


When Dg^Dy = 1, equation (C6) reduces to the expression obtained by Seide (ref. 52) for 
unstiffened shells. Note that equation (C6) gives the minimum value of p C y^ directly. 

The correlation curve in figure 18 was obtained by Seide (ref. 52). The open sym- 
bols represent results from the present investigation for unstiffened shells; values of 
Pcyl obtained from equations (C3) and (C6) were essentially the same. The agreement 
is quite good. Most of the scatter is for a g 45°; Singer (ref. 55) has shown that the cor- 
relation curve overestimates the ratio P/P C yi somewhat for this range of a. Results 
presented in reference 54 indicate the correlation curve can also be used to correlate 
data for orthotropic and ring-stiffened conical shells. Accordingly, results for ring- 
stiffened shells are also shown in figure 18 as solid symbols. Equivalent orthotropic 
stiffness properties were calculated from the formulas given in table I, which neglect 
eccentricity effects. 


The solid squares in figure 18 represent results from reference 56 (based on 
eq. (C6)) for Dgy/Dy 2 25 and agree quite well with the correlation curve. Results 
from the present investigation for Dp/Dy = 316 (ring details are given in fig. 12) are 
indicated by the solid circles for p C yj defined by equation (C6) and by the solid symbols 
with ticks for p C y^ defined by equation (C3). The scatter shown by these results sug- 
gest that the correlation curve becomes more approximate as Dg/Dy increases and that 
a family of curves would probably be required for proper correlation, particularly for 
larger values of the taper ratio ^1 - R^/R2 > 0.6^. Use of the more exact definition of 
Pcyl ( e 1* (C3)) resulted in larger values of p/p C yi and somewhat better correlation than 
the results obtained by the use of the approximate expression (eq. (C6)). 

Effects of ring or s tringer stiffening .- The relative magnitudes of the effects of ring 
and stringer stiffening are indicated in figure 19 in terms of the ratio of the buckling load 
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of the stiffened shell to that of the unstiffened shell and L/Rq. Flutter results for the 
same shells are given in figure 12. The open symbols are for ring stiffeners and the 
solid symbols are for stringer stiffeners; stiffener details are given in figure 12. The 
equivalent orthotropic stiffnesses were based on the stringer spacing at the small end of 
the cone. Stiffener size and spacing was such that the ring- or stringer -stiffened shells 
were essentially of equal weight. The vertical lines connecting the open symbols in fig- 
ure 19 represent the variation in the ratio of the buckling loads with cone angle a ; in 
general, the ratio decreased with increasing cone angle. The variations with cone angle 
for stringer stiffening was not significant. The advantage of ring stiffening (or orthotropic 
stiffening for D^/Dy > l) for hydrostatic pressure loading is obvious. The curves shown 
in figure 19 were obtained from the expressions for equivalent cylinders (eqs. (C3) and 
(C6)). The variation with a (for fixed l/Rj) was not significant; thus, only a single 
curve (for each equation) is shown. The qualitative effect of stiffening indicated by the 
equivalent cylinder equations is in reasonable agreement with the actual, calculated cone 
results. 


Elastically Supported Shells 

The results of the preceding sections show that from a theoretical standpoint the 
buckling of a simply supported unstiffened conical shell subjected to hydrostatic pressure 
is adequately described by equation (C6) in conjunction with the correlation curve shown 
in figure 18. Furthermore, the buckling of ring-stiffened shells is fairly well described 
by equation (C3) or (C6), and the results obtained, if not sufficiently accurate for final 
buckling criteria, should at least prove useful for preliminary calculations. On the other 
hand, the considerable scatter in available experimental data is such that experimental 
values range from 60 percent to 140 percent of the values predicted for simply supported 
edges (ref. 3); the scatter is attributed to variation in cone edge conditions and to initial 
imperfections. In the next three sections, the effects of variations in edge conditions are 
shown to have a significant effect on the buckling of conical shells. 

Effects of ax ial and rotational restraint.- Some effects of axial restraint (kj) are 
shown in figure 20, which gives the variation of the buckling pressure coefficient Cp 
with the ratio k^/E. Fairly large values of axial restraint (ki/E > O.OOl) are required 
to noticeably affect Cp, but as the restraint approaches complete constraint the effect is 
appreciable. The value of Cp for complete constraint is about 40 percent greater than 
the value for no axial restraint. The symbols shown in figure 20 represent buckling 
results for no axial restraint and complete axial restraint obtained previously by Thurston 
(ref. 24) who used a very accurate prebuckling state. The results of the present investi- 
gation are in excellent agreement with the results of reference 24. (See fig. 20.) This 
agreement suggests that the extra computation involved in considering the accurate pre- 
buckling state is not necessary, at least for hydrostatic pressure loading and the edge 
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conditions considered. Similar results and conclusions for no axial restraint and com- 
plete axial constraint were presented in reference 57. 

The variation of the effects of complete axial restraint with L /R 4 are shown in 

PuT7"r> = 0 

figure 21. The variation of the buckling ratio — - — increases rapidly with L/R* 

p v B =o ' 1 

up to about 2 and apparently would reach a limiting value of between 1.4 and 1.5 for very 
large values of L/Rj. More complete results on the variation of the effects of no axial 
restraint and complete axial constraint with L/Rj, Rj/h, and a are given in refer- 
ence 57 where similar trends are indicated. 

Figure 22 shows some effects of rotational restraint in terms of Cp as a function 
Rik 4 

of . The lower curve is for no axial restraint, and the upper curve is for complete 

D y 

axial constraint. The results show that increasing rotational restraint increases the 
buckling pressure about 8 percent for no axial restraint but only about 2 percent for com- 
plete axial constraint. This tends to confirm the results of other investigations (for 
example, ref. 11) that the effect of a given edge restraint depends on the magnitudes and 
types of other restraints acting on the shell edges. For complete axial constraint and 

^^4 greater than about 500, the restraint approaches the fully clamped condition. The 

square represents buckling results for the fully clamped condition previously obtained by 
Singer (ref. 58); the computation is in excellent agreement with the results of the present 
investigation. The results shown in figure 22 indicate that the large increase in buckling 
pressure for clamped edges (square) compared with the results for cylindrical simple 
support (circle) is primarily caused by axial constraint. The effect of axial restraint is 
sufficiently large to preclude neglecting it in comparisons of theoretical and experimental 
results. Indeed, improved predictions of experimental results have been shown to occur 
when the axial constraint present in the experiment is considered in the analytical com- 
putations (ref. 57). 

Effects of^circumferential and r adial restraint .- Some effects of finite circumferen- 
tial (k2) and radial ^3) elastic restraint are shown in figure 23, which gives the variation 
of Cp with k 3 /E for a sandwich shell subjected to lateral pressure and axial load such 
that the large end of the cone is free of stress. Shell structural properties are given in 
table II. The upper curve in figure 23 is for k 3 /E = 1 and the lower curve is for 
k 3 /E = kg/E; for both curves k 4 = k 4 = 0. The circle and square at k 2 /E = 1 show the 
slight effect of changing the circumferential restraint from complete constraint ^Vg = 0) 
to no restraint (Ny^g = 0) for no radial displacement. The two symbols (diamond and 
triangle) at kg/E = 0 show the effect of changing the circumferential boundary condition 
from Vg = 0 to Nygg = 0 for no radial restraint; the effect is in marked contrast to 
the effect when there is complete radial constraint. The reason for the large effect of 
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circumferential restraint as the radial restraint becomes weak is that the constraint 
v B = 0 prevents the inextensional form of buckling (rj = 2), which is the critical mode of 
buckling when all edge restraints are relatively weak (refs. 25 and 59). Since the critical 
value of rj becomes small as the edge restraint becomes relatively weak, the present 
Donnell type theory becomes more approximate. Hence it might be desirable to use more 
exact shell theories in analyses of buckling of shells with weak edge restraint. 

The sandwich shell considered in figure 23 has previously been analyzed by Cohen 
(ref. 25) who utilized Novozhilov shell theory (ref. 60) and very accurate prebuckling 
stresses and deformations. Although the main purpose of the analysis was to consider 
the effects of the edge restraint imposed by rings, results for simply supported edges 
were also presented and are indicated in figure 23 by the solid symbol. The result of 
the present investigation for simply supported edges is in good agreement with the result 
from the more rigorous analysis of Cohen, the difference being only about 6 percent. 

Effects of re straint of e nd rings .- Results for the sandwich shell considered in the 
previous section and a similarly loaded ring- stiffened shell are presented in figure 24 in 
terms of Cp as a function of the ratio of the radius of the tubular end ring (of circular 
cross section) to the radius of the shell at the large end (r/R) 2 - The ratio (r/R)i for 
both shells was constant and equal to 0.0312; structural details of the shells are given in 
table II. For the end rings, rj = 1.25 in. (3.18 cm), tj = T 2 = 0.125 in. (0.318 cm), and 
r 2 is varied. When the edge restraint is weak (small values of (r/R^), the shell tends 
to buckle inextensionally into two circumferential waves, but when the edge restraint 
becomes large, the shell buckles extensionally into seven circumferential waves. The 
intersection of the curves for r\ = 2 and rj = 7 gives the minimum base-ring radius 
ratio (r/R )2 required to suppress inextensional buckling; further stiffening would not 
have a significant effect on the buckling pressure. The minimum value of (r/R )2 to 
suppress inextensional buckling obtained from the present analysis is conservative in that 
it is approximately 35 percent larger than the value predicted by the more rigorous anal- 
ysis of reference 25 for the sandwich shell and 18 percent larger for the ring-stiffened 
shell. These differences are due primarily to use of Donnell theory and the membrane 
stress state in the present analysis and, also, to neglect of ring eccentricity. 

The results of figure 24 indicate that the present analysis gives fair predictions 
(compared with the more rigorous analysis of ref. 25) of the buckling pressure for ring- 
stiffened pressure -loaded shells. However, for axial load or stringer stiffening the 
neglect of eccentricity effects and the use of the membrane prebuckling state can lead to 
extremely inaccurate results. For example, in reference 61, which utilizes a Donnell 
type buckling analysis, it is shown that either eccentricity effects or the use of a very 
accurate prebuckling state can cause over a factor of 1.5 change in the buckling load. 
Results of the present analysis differed by about a factor of 2 from the results of 
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reference 61 for a clamped shell with a = 5° and external stringers. Thus, the present 
analysis should be used with caution for configuration and loading combinations sensitive 
to prebuckling and eccentricity effects. 
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VIBRATION RESULTS 
Simply Supported Shells 

Effects of conicity .- Some effects of coni city on the vibration characteristics of 
shells with cylindrical simple support are shown in figure 25. The decrease of (c u/a> r ) 2 
with increasing a (fig. 25(a)) and l/Rj (fig. 25(b)) is similar to the trends for 
buckling (fig. 17). Extrapolated results for zero cone angle are seen to be in excellent 
agreement with the cylinder results obtained from equation (43) of reference 62, indi- 
cated by the circles in figure 25. 

The minimum frequencies for isotropic conical shells have been correlated (ref. 63) 

Rl 

by use of the taper ratio 1 and the concept of equivalent cylinders. Results from 

r 2 

the present analysis for 5° = a = 60°, 1.26 = L/Rj ~ 6.28, and 107 = R i/h = 341 are 

compared with the correlation curve of reference 63 in figure 26 which gives the varia- 

Ri 

tion of w min/( w cyl) m i n wR h 1 “ An expression for the natural frequency of an 

equivalent cylinder applicable for unstiffened or ring stiffened shells can be 

obtained by substituting equations (C2) into the cylindrical vibration equation (eq. (43) of 
ref. 62). The result for the lowest frequency (m = 1) for a given value of 77 is 
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where js is defined by equation (C4). It should be noted that equation (43) of refer- 
ence 62 contains expressions multiplied by Zr, which can be either positive or nega- 
tive. These terms, which account for eccentricity (one-sided) effects, are omitted in 
equation (Dl). 

Equation (Dl) is a function of 77 , which must be varied to find ( w cyl) m i n * For 
unstiffened shells the expression for ( w cyl) m i n becomes (ref. 63) 
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As can be seen from figure 26, the results of the present analysis for unstiffened shells 
are in good agreement with the correlation curve. Calculations were also made for ring- 
stiffened shells for which Dg/Dy = 316 as indicated by the solid circles in figure 26. 

As can be seen, the equivalent cylinder vibration results for stiffened shells correlate 
better (fig. 26) than did the buckling results for similar shells (fig. 18). It should be 
noted that the equivalent cylinder concept has limited usefulness for frequencies 
u> > w m i n (ref. 63). As rj increases, the mode shapes become more confined to the 
region near the large end of the cone and, therefore, can deviate considerably from the 
sinusoidal shape of a vibrating simply supported cylinder. Thus, the difference in the 
frequencies for actual cones and equivalent cylinders increases as 77 increases; this 
trend was also noted by Weingarten (ref. 64) for ring-stiffened conical shells. 

Effects of applied loads .- Effects of hydrostatic pressure and axial load on the nat- 
ural frequencies of an isotropic conical shell with cylindrical simple support are shown 
in figure 27 in terms of (w/o > r )2 as a function of the ratio p/p cr or P/P C r f° r sev “ 
eral values of 77 . As can be seen from figure 27(a), (ayu> r )2 varies linearly with 
p/p cr a nd the slope of the curves increases with increasing 77. For the range of 77 
considered, which brackets the value for buckling (?) = 10 ), (w/a > r )2 decreases signifi- 
cantly with p/p cr . Weingarten (ref. 65) has presented results which also show signifi- 
cant effects of increasing p/p cr and, in addition, indicate that the frequency is nearly 
independent of p/p cr for 77 ^ 2 . 

Some effects of axial load are shown in figure 27(b) for the same shell considered 
in figure 27(a). In contrast to the results for hydrostatic loading, figure 27(b) reveals 
that (o>/w r )2 does not vary linearly with P/P cr ; this is attributed to the differences 
in the axial variation of Ny^ for hydrostatic loading and axial loading. For 
P/P cr < 0.5 the variation is nearly linear and the slope of the curves again increases 
with 77 . 
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Effects of rin g or stringe r stiffening .- The relative magnitudes of the effects of 
ring and stringer stiffeners on the minimum frequencies of conical shells are indicated 
in figure 28 in terms of the ratio of the frequency of the stiffened shell to the frequency 
of the unstiffened shell as a function of L/Rj. Stiffener size and spacing was such that 
the shells were of equal total weight; stiffener details and flutter results are given in fig- 
ure 12, and buckling results are given in figure 19. The equivalent orthotropic stiffnesses 
were calculated from the formulas given in table I and were based on the stringer spacing 
at the small end of the shell. The variation of the frequency ratio with a. was not sig- 
nificant for the range considered (5° = a g 60°); thus, an average value is presented. The 
solid curves represent results for actual conical shells, and the dashed curve represents 
results for equivalent cylinders; the results for equivalent cylinders are in excellent 
agreement with the results for actual conical shells. For the particular configuration 
studied, the effect of rings is to increase the minimum frequencies by a factor of 2 to 3. 
The addition of stringers has a slight effect, and may even lower the frequency. The 
effect on vibration results of adding stiffeners can be (and for the configuration studied 
is) considerably less than the effect for buckling (fig. 19). The reason for the difference 
is that for vibration the additional mass of the stiffeners can tend to offset the effects of 
increased stiffness. 


Elastically Supported Isotropic Shells 

The results of the preceding sections show that simply supported unstiffened or 
ring-stiffened conical shells are fairly well described by equation (Dl), or equation (D2) 
for unstiffened shells in conjunction with the correlation curve shown in figure 26. 
Agreement between theoretical and experimental results is generally good (for example, 
refs. 63 to 66); discrepancies that exist are attributed to effects of edge conditions 
(ref. 63). In the next two sections it is shown that the effects of variations in edge condi- 
tions can have a significant effect on the vibration characteristics of conical shells. 

Effects of axial and r otational restraint .- Some effects of variations in elastic axial 
restraint and rotational restraint (k^ are shown in figure 29. Figure 29(a) shows 
the variation of (oyo> r )2 with kf/E for no rotational restraint, and figure 29(b), shows 

9 Rik4 

the variation of (a>/cn r ) 4 with ■ ^ — for no axial restraint and complete axial constraint. 

As can be seen from figure 29, the effects of axial and rotational restraint on free vibra- 
tions are similar to the trends for buckling. (See figs. 20 and 22.) The effect of axial 
restraint appears to be sufficiently large to preclude neglecting it in comparisons of theo- 
retical and experimental results. For example, in reference 66 the agreement between 
theoretical and experimental frequencies was improved by relaxing the axial constraint 
condition in the theoretical calculations. On the other hand, increasing axial restraint 
improves the agreement between theoretical frequencies and the experimental frequencies 
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9 

presented in reference 65 as shown in figure 30 in terms of (co/o> r ) as a function of rf 
for P/P cr = -0.446, 0, and 0.446; p cr is the theoretical buckling pressure. The solid 
curves are for the conical simple support boundary conditions Nyg = Vg = Wg = Myg = 0 
and agree with the theoretical results presented in reference 65. The dashed curves are 
for the cylindrical simple support, Vg = Vg = Ugg = Myg = 0, and indicate slightly higher 
frequencies than the solid curves. The long-dash — short-dash curves are for the same 
support conditions as the dashed curves except that k^/E = 0.03; this amount of axial 
restraint was chosen so that the theoretical and experimental results would agree for 
^min = ® and p = 0. As can be seen from figure 30, inclusion of some finite axial 
restraint in the theoretical calculations improves the agreement between theoretical and 
experimental frequencies. 

Effects of circumferential and radial restraint .- Some effects of variations in elas- 
tic circumferential restraint (k 2 ) and radial restraint (k 3 ) are shown in figure 31 in terms 
of (a>/u> r )2 as a function of the ratio k 3 /E. The upper curves are for r) - 10 and the 
lower curves are for rj = 6 . The open circle and square at kg/E = 1 
slight effect of changing the circumferential restraint from complete constraint (Vg = 0 ) 
to no restraint (Nygg = 0) for no radial displacement. The results of the present inves- 
tigation are in good agreement with the results of reference 63 as indicated by the solid 
symbols in figure 31; the results for Ugg = 0 (present investigation) are slightly higher 
than the results for Wg = 0 (ref. 63). The trends exhibited by the curves for k 2 /E = 1 
(vg = 0 ) and k 2 ^E = k 2 /E as k 2 /E decreases are similar to the trends for buckling 
(fig. 23). Results for other values of rj would exhibit similar trends although the mag- 
nitude of the change in frequency due to changing the edge restraint would vary with r j. 

For example, in reference 63 it is shown that changing the edge restraint from Vg = 0 
to N y 0 g = 0 (for Wg = 0 ) changes the frequency by about 50 percent for rj = 2, whereas 
for r) § 10 the change was about 2 percent or less as indicated in figure 31. 


(u HB = 0j show the 
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TABLE I 


ELASTIC CONSTANTS FOR RING OR STRINGER 
STIFFENED CONICAL SHELL 
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TABLE II 


STRUCTURAL PROPERTIES OF CONICAL MODELS 
STUDIED IN REFERENCE 25 


Sandwich shell 

Face sheets 

Core 

E = 10.5 X 10 6 psi (72.5 GN/m2) 

E = 0 

G = 3.98 X 10 6 psi (27.5 GN/m2) 

G = 0 

/i = 0.32 

ju = 0 

h = 0.020 in. (0.051 cm) 

h = 0.50 in. (1.27 cm) 

Ring-stiffened shell 

Skin 

Ring stiffeners 

E = 6.5 X 106 psi (44.8 GN/m2) 

E r = 6.5 X 10 6 psi (44.8 GN/m2) 

G = 2.4 X 10 6 psi (16.6 GN/m2) 

G r = 2.4 X 10 6 psi (16.6 GN/m2) 

p = 0.35 

p R = 0.35 

h = 0.051 in. (1.30 cm) 

I R = 0.00414 in 4 (0.172 cm 4 ) 


J R = 2i R 

Ar = 0.0635 in^ (0.410 cm2) 
dR = 1.8 in. (4.57 cm) 
z cR = ^-164 in. (0.418 cm) 
z sR = ^-236 in- (0-600 cm) 



(b) Force and moment resultants. 



(a) Configuration; coordinates and displacements. 


(c) Loads. 


Figure 1.- Configuration and notation for coordinates and displacements, force and moment resultants, and loads. 
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Figure 2.- Variation of flutter parameter \ with wave number n for isotropic conical shell with cylindrical simple support. 

L/ R i = 1; Rj/h = 200; a = 20°; p = P = 0. 


57 




a, deg 


(a) Cone angle. 
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(b) Length-radius ratio. 

Figure 3.- Some effects of conicity on flutter of isotropic conical shells with cylindrical simple support. Rj/h = 200; 
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Figure 15.- Some effects of rotational restraint on flutter of isotropic conical shell with no axial restraint and complete axial constraint. 

L/Rj = 4 ; Rj/h = 200; a = 30°; p = P = 0; g a = 0; v B = u H g = 0. 


k 2 /E 


F I utter 


-k 2 /E = k 3 /E 


No flutter 


T) = 8 


□ v B =0 


■ V B=° 


0 ( U HB = C 


Q yB - 0 


Figure 16.- Some effects of elastic circumferential (k 2 ) and radial (k 3 ) edge restraint on flutter of isotropic conical shell 
V R 1 = 3 - 5 ; R l/h =107; a = 200; p = P = 0; g a = 0; V = M vB = 0. 
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Figure 18.- Correlation of results for buckling of equivalent cylindrical shells with conical shells with cylindrical simple support 

subjected to hydrostatic pressure. P = 0. 
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Figure 19.- Some effects of ring and stringer stiffening on buckling of conical shells with cylindrical simple support subjected to 
hydrostatic pressure. Rj/h = 341; P = 0; 5° ^ a ^ 60°; ring and stringer details are shown in fig. 12. 
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Figure 20.- Some effects of elastic axial restraint on buckling of isotropic truncated conical shell subjected to hydrostatic pressure 

L/Ri = 2; Rj/h = 500; a = 30°; P = 0; v B = u HB = M yB = 0. 





Figure 22 - Some effects of rotational restraint on buckling of isotropic truncated conical shell subjected to hydrostatic pressure. 

L/ R l = 1; R ly /h = 707- a = 45°; P = 0; v B = u HB = 0. 


Investigation 

Edge condition o present 



Figure 23.- Some effects of elastic circumferential (k 2 ) and radial ^ 3 ) edge restraint on buckling of truncated conical sandwich shell 
subjected to lateral pressure. L/Rj = 0.72; a = 60°; Vg = Myg = 0; shell details are given in table II. 
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(a} Cone angle. 



(b) Length- radius ratio. 

Figure 25.~ Some effects of conicity on natural frequencies of isotropic conical shells with cylindrical simple support. 

Rj/h = 341; n = 10; p = P = 0. 




Figure 26 - Correlation of minimum frequencies of equivalent cylindrical shells with conical shells with cylindrical simple support, p P 0. 



(a) Hydrostatic pressure. fW Axial load. 

Figure 27.- Some effects of applied loads on natural frequencies of isotropic conical shell with cylindrical simple support. 

L/Rj = 4; Rj/h = 200; a = 30°. 
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Figure 28.- Some effects of ring and stringer stiffening on minimum natural frequencies of conical shells with cylindrical simple support. 
Rl/h = 341; 5° = a ~ 60°; ring and stringer details are shown in fig. 12; p = P = 0. 



(a) Effects of axial restraint with no 
rotational restraint. 


(b) Effects of rotational restraint with no axial 
restraint and complete axial constraint. 


Figure 29.- Some effects of elastic axial restraint (k]J and rotational restraint (k^) on natural frequencies of isotropic conical shell. 
L/Ri = 4; Ri/h = 200; a = 30°; n = 10; v B = u HB = 0; p = P = 0. 
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(a) p/p cr = -0.446. 


(b) p/p cr = 0. 


(c) p/p cr = 0.446. 


Figure 30.- Comparison of theoretical and experimental vibration results for isotropic conical shell L/Ri =3 5- Ri/h = 107- 

a = 20°; P = 0; v B = M yB = 0. ’ ' 
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Figure 31.- Some effects of elastic circumferential (k 2 ) and radial (k 3 ) edge restraint on natural frequencies of isotropic conical shell. 

L/Rl = 3.5; Rj/h = 107; a = 20°; p = p = 0; V B = MyB = 0. 
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The lengthy expressions developed in the analysis and required for complete 
documentation are included in a "Supplement to NASA TN D-5759." 

Copies of this "Supplement to NASA TN D-5759" will be furnished upon 
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